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ABSTRACT 

Previous studies of the Na I D interstellar absorption line doublet have shown 
that galactic winds occur in most galaxies with high infrared luminosities. How- 
ever, in infrared-bright composite systems where a starburst coexists with an 
active galactic nucleus (AGN), it is unclear whether the starburst, the AGN, or 
both are driving the outflows. The present paper describes the results from a 
search for outflows in 35 infrared-faint Seyferts with 10 9 ' 9 < Lir/L < 10 11 , or, 
equivalently, star formation rates (SFR) of ~0.4 - 9 M yr _1 , to attempt to 
isolate the source of the outflow. We find that the outflow detection rates for the 
infrared-faint Seyfert Is (6%) and Seyfert 2s (18%) are lower than previously re- 
ported for infrared-luminous Seyfert Is (50%) and Seyfert 2s (45%). The outflow 
kinematics of infrared-faint and infrared-bright Seyfert 2 galaxies resemble those 
of starburst galaxies, while the outflow velocities in Seyfert 1 galaxies are signifi- 
cantly larger. Taken together, these results suggest that the AGN does not play a 
significant role in driving the outflows in most infrared-faint and infrared-bright 
systems, except the high-velocity outflows seen in Seyfert 1 galaxies. Another 
striking result of this study is the high rate of detection of inflows in infrared- 
faint galaxies (39% of Seyfert Is, 35% of Seyfert 2s), significantly larger than in 
infrared-luminous Seyferts (15%). This inflow may be contributing to the feed- 
ing of the AGN in these galaxies, and potentially provides more than enough 
material to power the observed nuclear activity over typical AGN lifetimes. 

Subject headings: galaxies: active - - galaxies: Seyfert - - galaxies: quasars: 
absorption lines — line: profiles — ISM: jets and outflows — ISM: kinematics 
and dynamics 
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INTRODUCTION 



Galactic-scale gas outflows appear to play an important role in the evoluti on of the uni- 
verse , possibly providing an explanation for a number of cosmological questions (jVeilleux et al 
20051 and references therein). These galactic outflows are frequently observed in local galax- 
ies with high star formation rates and/or an active nucleus, often extending on a scale of 
order of a few kiloparsecs or larger. 

The high frequency of outflows measured in high-redshift, actively star- forming ob- 
jects such as Lyman break galaxies suggests that ou tflows are a common stage in galactic 
evolu tion, with wind velocities decreasing over time (jShapley et al.ll2003l : iFerrara fc Ricotti 
20061 ). Such outflows may be responsible for the deficit of baryons seen in our own Milky 
Way galaxy and the mass-metallicity relation in external gal axies, if winds are capable of 



preferentially electing metals into the intergalactic medium (jLarsonl Il974l ; iGarnettl 12002 



Tremonti et al.ll2004l ). The intergalactic medium can be heated as well as enriched with 
metals by these galactic outflows, since energy inputs into th e IGM on the order of 10 56 
ergs have been measured from galaxies with strong outflows (ICroston et aJj 120081 ) . Addi- 
tionally, outflows may quench star formation b y heating up cold gas and ejecting it from 
the host (jBinneyl 12004 ; IScannapieco et al.ll2005l ). Specifically, AGN-powered outflows have 
been proposed as the cause of the drop in AGN lu minosity at low redshift and the cutoff at 



high l uminosity of the gal axy luminosity function (jSomerville fc Primack! 1 19991 ; ICole et al. 



2000l ; iThacker et al.ll2006l ). These winds may also limit black hole and spheroid growth, 



and be responsibl e for the observed tight correlation between black hole mass and galac - 
tic spheroid mass (IFerrarese &; Merrittll2000l ; iMcLure &; Dunlopl 120021 ; iGranato et al.l I2004J ) . 
Finally, winds driven by AGN and/or starbursts may help remove enough nuclear material 
with low angular momentum early on in the evolution of gas-rich systems to aid in the forma- 
tion of large disks with high specific angular mom entum, more in line with the observations 
( IBinnev et al.ll200ll ; iKormendv fc Kennicuttl 12004 ). 



Studies performed on actively star-forming galaxies, near and far, have revealed that 



outflow energeti cs increase wi t 
formation rate (IMartinl 12005 



hereafter RVS05b; ISato et al. 



i infrared (8-1000 a m) luminosity, L IR , or equivalently, star 



; iRup 
2009; 



Weiner et al. 



ke et all l2005aL he reafter RVS05a; iRupke et all [2005b 



2009i ). Outflows have also been detected in 



comp osite galaxies with both starbursts and AGN (e.g.. llrwin &; Saikiall2003l ; iHota fc Saikia 
20061 ). Seyfert 1 composites sho w significantly high er outflow velocities than pure star- 
bursts and Seyfert 2 composites ( Rupke et al. 12005c, hereafter RVS05c). Hydrodyna mical 
simulations of starburst-driven outflows (e.g.. IThacker et al.ll2006l ; ICooper et al.ll2008l ) have 
difficulties reproducing the very high velocities detected in Seyfert 1 composites. The AGN 
in these objects thus appear to play a significant role in driving these outflows, at least on 
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small scales (ICrenshaw et al.ll2003bl ). The situation in Seyfert 2 composites is more ambigu- 



ous - the AGN or starburst or both could be powering the outflow in these systems ( jSchiano 



19851 ). 



The focus of this paper is on "pure" Seyfert galaxies with weak infrared starbursts. Our 
sample consists of 35 galaxies that are faint in the infrared (10 9 ' 9 < Lir/Lq < 10 112 ), equally 
split between Seyfert Is and Seyfert 2s. We compare the results from our study with those 
from previous studies of starburst and Seyfert ultraluminous infrared galaxies (ULIRGs; 
L IR /L > 10 12 ) and luminous infrared galaxies (LIRGs; 10 11 < L IR /L < 10 12 ) to isolate 
the role of the AGN in powering these outflows. As in previous studies, the IR-faint objects 
were observed in the Na I D AA5890,5896 doublet absorption feature. This Na feature has a 
low ionization potential (5.1 eV), so it probes neutral gas, and it can be used to study the 
ISM due to its high interstellar abundance. All objects in the sample have z < 0.05, and 
thus the Na I D feature is found in the optical. Velocity components in Na I D absorption 
are unambiguous indicators of outflowing (blueshifted velocities) or inflowing (redshifted 
velocities) g there must be a continuum source behind the absorber. 

The organization of this paper is as follows. The sample is described in Section [2J 
including methods of selection and sample properties in terms of redshifts, star formation 
rates, and spectral types. The observations are discussed in Section [3j Section 0] describes 
the line fitting, the derivations of the velocities and column densities, and estimates of the 
expected stellar contributions to the measured Na I D absorption. In Section \5\ we present 
the results on outflows, starting with the Seyfert 2s and followed by the Seyfert Is. Next, 
Section [6] describes the results on inflows in the same fashion as Section [5j This is followed 
by a discussion of the implications of our findings in Section [7J including comparisons with 
previous studies, dynamical estimates, and possible connections between inflows and nuclear 
structures. We conclude in Section [S] with a summary of our findings. 

All calculations in this study assume present-day cosmological parameter values, with 
H = 75 km s _1 Mpc" 1 , Q m = 0.3, and = 0.7. All wavelengths quoted are vacuum 
wavelengths. 



2. SAMPLE 

The objects in the sample were selected using three main criteria: They had to be 
optically classified as Seyferts and had to be faint in the infrared (Lir < 10 11 L & ). All 
objects in the sample were also selected to have z < 0.05 to focus on the brightest sources 
and obtain reliable measurements. 
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Most galaxies in the sample are well-studied Seyfert galaxies with extensive ancillary 



data (e.g. IWhittldll992al ; iNelson &: Whittle! Il995l ). Roughly equal numbers of Seyfert 1 and 



Seyfert 2 galaxies were chosen to allow meaningful comparisons with the study of RVS05c 
on ULIRG Seyferts. The bolometric luminosities of the AGN in the present sample are 
well matched to those in the composite systems of RVS05c. This allows us to make direct 
comparisons between the two samples. Basic properties for the galaxies in our sample can be 
found in Table HJ The last column in that table lists the references used for the compilation. 



2.1. Redshifts 



Not all objects in this survey have well-constrained publishe d redshift measurements 



so preference was given to values from stellar absorption lines (INelson fc Whittle 



1995 



Wegner et al.ll2003 ) when p ossible, followed by redshifts taken from HI lines ( INelson fc Whittle 
)95l ; ISpringob et al.ll2005l ). Redshifts measured from data in the present survey were used if 
neither of the preceding were available, but were given low priority due to systematic errors 
(~ 10~ 4 ) associated with these measurements (e.g. it is possible that part of the narrow line 
region (NLR) is outflowing so the velocities based on the narrow emission lines are not neces- 
sarily representative of the systemic velocity). In this case, IRAF was used to measure central 
wavelength of various emission lines, including [O I] AA6300, 6364, [N II] AA6548, 6583, and 
[S II] AA6716, 6731, and redshifts were calculated from those values. 



2.2. Star Formation Rates 



To estimate the star formation rates in these objects, we use their far-infrared luminosity, 
L(40- 120 Atm), under the as sumption that the AGN does not contribute significantly to 
LpiR (jSchweitzer et al.ll2006l ). This is pr obably a good first -order approximation although 
AGN may contribute ~10-20% to Lfir (INetzer et al.l 120071 ). Lfir is calculated using the 
prescription of Sanders & Mirabel (1996), using data from the IRAS Faint Source and IRAS 
Point Source Catalogs. Lfir is then used to calculate SFR using the relation in Kennicutt 
(1998), SFR = Lpm/ (5.8xlO 9 L ) M yr _1 . The far-infrared luminosities and star formation 
rates of the objects in our sample are listed in Table [TJ 
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2.3. Spectral Types 

Spectral classifications were taken from the NASA Extragalactic Database^ and con- 
firmed through examination of emission lines in the spectra and references in the literature. 
Seyfert 1-1.5 were labeled Seyfert Is for the purpose of this study, while Seyfert 1.8 - 2 were 
labeled Seyfert 2s. This is in keeping with the classification scheme of RVS05c and allows for 
comparisons of our results on Seyfert 1 and Seyfert 2 galaxies with those of previous studies. 

3. OBSERVATIONS 

All observations were taken over the course of three different observing runs on the Kitt 
Peak 4-meter telescope. Observing runs, exposure times, and slit position angles for all ob- 
jects are listed in Table [TJ All data were taken using the Ritchey- Chretien Spectrograph with 
a moderate-resolution grating, KPC-18C, in the first order along with the GG-475 blocking 
filter. The wavelength range was 1700 A, allowing measurements of both the Na I D absorp- 
tion doublet and Ha + [N II] emission complex in the same exposure. This combination 
provides an average resolution of 85 km s _1 with a l'/25 slit. Median signal-to- noise ratio 
near Na I D was 85 per A and ranged from 23 to 337 per A, with seeing ranging from ~1 to 
2" on average. 

4. ANALYSIS 
4.1. Line Fitting 

Extraction and fitting procedures were performed as described in the original study 
(RVS05a), but the basic procedure will be outlined here; the method uses assumptions 
similar to those in a curve of growth analysis but does explicitly fit the line profiles. The 
spectral extraction was done while keeping a constant physical aperture size of ~3 kpc; this 
aperture size was selected to match the seeing and therefore reduce host galaxy contamination 
as much as possible. The reduction and calibration were performed using standard IRAF 
reduction packages. HeNeAr lamps were used for wavelength calibration and stars such as 
G191-B2B and BD+25 4655 were used as flux standards. 

Once the spectra were reduced and calibrated, the region containing the Na I D absorp- 
tion doublet, as well as the neighboring He I A5876 emission line, was isolated. A Levenberg- 



http : //nedwww. ipac . caltech . edu/ index .html 
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Marquardt fitting routine was used to fit one to two components to the He I emission and 
the Na I D absorption doublet line profile. For emission lines, the formula used to fit the 
intensity of the line was: 

/(A) = 1 + Ae« x - x ^( Xob » 2 (1) 

where the amplitude, A, the Doppler parameter, b (b = V2a), and the central wavelength, 
Ao, were free parameters. For absorption lines, the formula used to fit was: 

/(A) = 1 - C f (X) + CV(A)e~ Tl(A) - T2(A) , (2) 

where Cf is the covering fraction of the line (see Section l7.1.2p . and T\ and r 2 represent the 
optical depths of each of the two lines of the Na I D doublet; T\ corresponds to the line at 
5896 A, and r 2 = 2t\. In general, the optical depth is given by: 

Ti {\) = r . i0e -(A-A ! ,o) 2 /(A,ofe/ C ) 2 . ( 3 ) 

Equations 2 and 3 were used to fit each line of the Na I D doublet. For each Na I D doublet 
pair, the covering fraction, Cf, central optical depth and wavelength of the red line {t\ and 
A ,i), and line width, b, were free parameters in the fit. Covering fraction was assumed to 
be independent of wavelength (and the same for both doublet lines), and the optical depths 
of the two doublet lines are fixed by atomic physics (as stated above: r 2 = 2n). Because 
of these constraints, and the assumption of a Gaussian in optical depth space, both Cf and 
r can be constrained simultaneously for each doublet. It should be noted, however, that 
there is still some anticorrelation between Cf and r in the fit due to uncertainty in line 
shape and relative intensities. In the case of multiple, Gaussian velocity components, each 
velocity component arises in both lines of the doublet, so that Cf and r can be constrained 
independently for each velocity component. 

Initial guesses for these parameters were supplied to the fitting program and a best fit 
was determined and plotted. Either one or two components were used to fit the He I emission 
line. Unless an absorption fit with two separate velocity components (corresponding to two 
separate Na I D doublets) was considerably better than a fit with one velocity component, 
only one velocity component (one Na I D doublet) was used to fit the absorption features. 
Only six objects showed convincing evidence for a two velocity components in absorption. 
Examples of the deblending of the Na I D feature for two Seyfert 1 and two Seyfert 2 objects 
can be seen in Figure [TJ After fitting was performed, a Monte Carlo simulation with 1000 
iterations was then run, in order to determine the parameter errors. 

Particular care is needed when fitting the He I emission at 5876 A, since the He I emission 
line in these objects is generally quite broad and has multiple components. So, for those 
objects showing He I emission lines, the Ha + [N II] complex was also examined to verify 
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the accuracy of the fit to He I . An interpolation script was written in IDL to subtract out 
the [N II] AA6548, 6583 features and thus isolate the Ha line. The resulting Ha line was 
generally fit using two components, one narrow and one broad. The He I + Na I D complex 
was then refit, but this time parameters for He I emission lines were constrained to be 
those of the Ha line. The results from the fits with fully fixed parameters ("fixed") were 
then compared to those from the fits with free-floating parameters ("free"); agreement is 
observed within 2-5%. Again, the errors on these parameters were determined from the 
Monte Carlo simulations after all fits were completed. The reduced x 2 f° r the two types of 
fits are compared in Figure [21 There are fewer points farther away from a reduced \ 2 °f 
unity for the free-floating fit. We therefore choose to use the results from the free-floating 
fits in the following discussion (this is further justified in the next section). 



4.2. Velocities 

Once fitting was done, velocities for these Na I D absorption lines were calculated, 
relative to the systemic velocities using redshifts as outlined in Section 12.11 In the following, 
Av is the velocity difference between the Na I D velocity centroid measured in the fit and 
the systemic Na I D velocity. Errors were propagated through both from the fit and the 
redshift measurement. A "maximum" velocity of Av max = Av - FWHM/2, which takes 
into account both the shift and width of the absorption line, was also computed. Values 
of Av and Av max were compared for values resulting from the free He I fit and the fixed 
He I fit in order to determine the impact of fixing the He I parameters to those of Ha. Plots 
of those comparisons are displayed in Figure [31 The agreement is good, with only three 
objects showing values that are not equal within the uncertainties. Of those, only one has 
l^^frcc — Affixed I > 50 km s _1 , which would change its outflow or inflow classification. In 
all three of those outlying cases, the He I line appears to be broader than the Ha line. This 
is plausible, since all of these objects are Seyfert 1 galaxies and so collisional and radiative 
transfer effects in the broad line region may affect the He I profile in a different way than the 
Ha profile. As there is little significant difference between the free-floating and fully-fixed 
fits, our adoption of the free-floating fits in Section 14.11 seems appropriate. 

After velocities were calculated, a detection criterion of Av < —50 km s _1 was used to 
determine whether or not an object possesses an outflow. Similarly, Av > 50 km s _1 was 
used as the cutoff for the detection of an inflow. A 2-cr threshold in measurement uncertainty 
(\Av\ > 28 (Av), where 8(Av) is the measurement error) was also required for inflow/outflow 
classification, but this did not eliminate any objects from being classified as showing inflow 
or outflow. As outlined in the initial study (RVS05a), the 50 km s _1 cutoff accounts for 
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potential contamination from errors in wavelength calibration (<10 km s _1 ), fitting (~10 
km s -1 ), redshift measurements (~ 10-15 km s" 1 ), and the possibility of small blue- and 
redshifts due to gas in rotation. Fit parameters and derived velocities are listed in Table [2J 



4.3. Column Densities 

The optical depth r and Doppler parameter b as determined by the fitting function were 
used to calculate column densities of Na I D and H along the line of sight. As described in 
Spitzer (1978), the column density of Na I D in cm -2 is given by: 

iV(Na I D) = , (4) 

V ; 1.497 x lO" 15 Ai/i' V ; 

where fi is the oscillator strength, Ai is the rest frame vacuum wavelength (A), and t± jC is 
the central optical depth of the Na I D A5896 line. The Doppler parameter, b, is in units of 
km s _1 . Values of fi = 0.3180 and Ai = 5897.55 A were taken from Morton (1991). 

To properly calculate the column density of hydrogen, we must correct for the effects of 
dust depletion and ionization. This process is also outlined in RVS05a, and uses empirical 
results that assume Galactic depletion and an ionization fraction of 0.9. The formula used 
is: 

N(R) = iV(Na I)(l - y)- l 10- {a+b) (5) 

Column densities are again given in cm -2 , y is the ionization fraction, and a and b are the Na 
abundance and depletion onto dust in the object for which calculations are being performed, 
respectively. Whereas a near-IR luminosity-metallicity relation was used to determine Na 
abundance in the previous RVS05b and RVS05c studies, we here assume that these objects 
have solar Na abundance. Calculated column densities for each object can be found in 
Table d 



4.4. Stellar Na I D Contribution 

Stellar Na I D absorption may contribute to the observed Na I D feature. Stellar 
features from other elements may also contaminate the measurements if they are located 
close to the Na I D feature. The original study used a scaling relation between Na I D and 
Mg I b equivalent width to determine the stellar Na I D contribution, as Na and Mg are 
created in a similar fashion (RVS05a). However, the spectral range for our data begins at 
roughly 5500 A and so does not include that Mg I b triplet, so we must use the other method 
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outlined in that study: stellar population synthesis models. As in RVS05a, we used the 
population synthesis code (Sed@.0 @) of Gonzalez Delgado et al. (2005), which combines a 
young (40 Myr), instantaneous burst stellar population with an old (10 Gyr) population. One 
model uses a stellar mass ratio of 10%/90% for young versus old populations, and the other 
uses l%/99%. In order to enhance weak stellar features, the spectra were boxcar smoothed 
by 150 km s _1 and convolved with a a = 200 km s" 1 Gaussian. Our spectra were emission 
subtracted, using fitting parameters to remove He I and leave only the Na I D absorption. 
They were then overlaid with the stellar population synthesis models. 

Almost all objects with suspected outflows or inflows showed absorption that was much 
deeper than the stellar absorption alone. The extreme examples of Mrk 352 and Mrk 6 are 
shown in Figure HI In Mrk 352, the results from stellar population synthesis show a high 
likelihood that the bulk of the Na I D absorption is stellar. In contrast, Mrk 6 shows a 
definite high-velocity outflow. 



5. RESULTS: OUTFLOWS 
5.1. Seyfert 2 

Figure [5] shows the He I and Na I D complexes for all of the 17 Seyfert 2 galaxies in 
this study, plotted on a velocity scale based on Na I D A5890 systemic. Tables [3] and H] list 
the measured properties for each Na I D outflow velocity component and overall outflow 
averages, combined with data for Seyfert ULIRGs from RVS05c. 

Out of a total of 17 Seyfert 2 galaxies, we detected blueshifted Na I D absorption with 
Av < —50 km s _1 in 3 objects, for an outflow rate of 18 ± 9%. No objects were ruled out by 
the 2-<t criterion in measurement uncertainty, which requires \Av\ > 25 (Av) for a detection 
(see Section 14721 . Of those objects which fit the outflow criteria, only one, NGC 7319, showed 
convincing evidence for two velocity components in absorption. The implications of this low 
detection rate are discussed in Section ITU 

Velocities for these objects were all within 20 km s~ x of each other, ranging from —130 
to —148 km s -1 . The maximum velocity, Av max = Av — FWHM/2, showed a much larger 
range (due to the broader range of Doppler parameters), from —176 to —504 km s -1 , with 
NGC 7319 showing the largest Av max (Table Hj). We further discuss these results in Section 

o 



2 http : //www. iaa. es/~rosa/ and http://www.iaa.es/" mcs/sed@ 
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Previous studies have shown that it is common for Seyfert 2 galaxies to have asymmetric 
emission lines with extended blue wings. These asymmetries are generally believed to be 
the result of outflowing gas in the narrow line region of the AGN, and outfl owing gas in 



the NLR has been spatia l ly resolved in both Seyfert 1 and Sey fert 2 galaxies (lWalkeiill968 



Wilson & Ulvestadl Il987i lEbstein et all Il989l : IVeilleuxi Il991al ffl iRuiz et all l2005h . This 



emission line asymmetry, however, is not an unambiguous indicator of outflow in the way 
absorption is since one cannot differentiate between outflowing line-emitting material in 
front of the nucleus and infalling material on the far side of the nucleus. Blue-wing emission 
asymmetries can only be interpr eted as outflow if we are abso lutely sure we are not seeing 
the opposite side of the object ([Rodriguez- Ardila et al.ll2006l ). The primary focus of our 
paper is blueshifted Na I D absorption, but it is of interest to look for blue-wing emission 
line asymmetry in our Seyfert 2 galaxies to see if it is at all correlated with our findings on 
Na I D outflows. In RVS05c, 75% of Seyfert 2 nuclei showed blue emission-line asymmetry 
(BELA), compared to a 45% Na I D outflow detection rate. 

We have used [N II] AA6548, 6583 to look for BELA in our objects (the high-ionization 
[O III] A5007 line most commonly used in BELA detection is not within our wavelength 
range). The [N II] AA6548, 6583 lines were separated from Ha by taking the blue wing of 
the [N II] A6548 line and the red wing of the [N II] A6583 line, and plotting them together 
about Av = 0. The two wings were then scaled appropriately to match each other in intensity, 
and BELA determination was done by visual inspection. 5 out of 17 Seyfert 2 galaxies (29%) 
showed evidence for BELA. The most obvious case of BELA is NGC 5252. Of those which 
show BELA, only one (NGC 7319) has been determined to have an Na I D outflow, and one 
(Mrk 348) actually shows a Na I D inflow (Section 16. ID . though it is the least asymmetric of all 
the BELA objects. Those five objects, along with a comparison object with a Na I D outflow 
but no BELA (Mrk 622), are shown in Figure[6j Another object (NGC 3786; not shown in the 
figure) shows prominent red-wing asymmetry; this object presents an Na I D outflow rather 
than an inflow. Overall, we find no obvious correlation between emission and absorption 
signatures of outflow, although the statistics are poor. Moreover, we cannot formally rule 
out the possibility that the difference in outflow and BELA detection rates is simply due to 
our use of [N II] as a BELA probe since BELAs a re more often s een in high-ionization lines 
like [O III] than in low- ionization lines like [N II] (I Veilleuxl 1 1 9 9 1 al ) . 



5.2. Seyfert 1 



Figure [7] shows the He I and Na I D complexes for all of the 18 Seyfert 1 galaxies in 
this study, plotted on a velocity scale relative to Na I D A5890 systemic. As for the Seyfert 
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2 galaxies, Tables [3] and H] list the measured properties for each Na I D outflow velocity 
component and overall outflow averages. 

Only 1 of the 18 Seyfert 1 galaxies was found to have blueshifted absorption with Av 
< —50 km s -1 , and as with Seyfert 2 inflows, the 2-a criterion did not eliminate any objects 
(Section 14.21) . The outflow detection rate is therefore only 6 ± 6%. 

Mrk 6 was the only Seyfert 1 galaxy determined to have an outflow. Two velocity 
components were detected, each with a velocity higher than those for the Seyfert 2 galaxies: 
one component with —229 km s _1 and one with —1024 km s _1 , giving a Av max = —1037 
km s _1 (Table H]). See Section I7TTJ for a discussion of these results. 

6. RESULTS: INFLOWS 

While only a few objects in our sample showed blueshifted Na I D absorption, a large 
percentage of objects in each Seyfert group showed redshifted Na I D absorption. We consider 
this redshifted absorption to be an unambiguous indicator of inflow, in the same way that 
blueshifted absorption is an unambiguous outflow indicator, since the continuum source (in 
this case the galaxy nucleus) must be behind the absorber. This section describes the results 
on inflow detection rates and kinematics. Table lists the measured properties of these 
inflows in individual objects. The implications of these results are discussed in Section 17.21 

6.1. Seyfert 2 

The same general detection criteria used for determining outflows were used for deter- 
mining inflows: We required a redshifted absorption of Av > 50 km s _1 , along with the 
2-a criterion (Section 14.21) . Out of the 17 Seyfert 2 galaxies, 6 showed Na I D inflow, for a 
detection rate of 35 ± 11%. 

Velocities for these six objects spanned the range from 51 to 127 km s _1 . The maximum 
velocity, where in the case of inflows Av max = Av + FWHM/2, ranged from 155 to 352 
km s -1 . The largest Av was measured for UGC 3995, although NGC 3362 has a slightly 
larger Av max than UGC 3995 due to the difference in Doppler parameter. The average 
properties of these inflows are listed in Table [6] and discussed further in Section 17.21 
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6.2. Seyfert 1 

Of the 18 Seyfert 1 galaxies, 7 showed Na I D inflow, for a detection rate of 39 ± 12%, 
i.e. similar to that of the Seyfert 2s. Again, no objects were ruled out by the 2a criterion 
(Section [472]). 

The Seyfert 1 inflow velocities spanned a similar range to those of the Seyfert 2 velocities, 
from 67 to 138 km s -1 . The Seyfert 1 Av max values showed a broader range than that of 
Seyfert 2, from 169 to 507 km s -1 , and while Mrk 1126 has the highest Av, NGC 7603 has 
the highest Av max . Again, the average inflow properties are listed in Tableland the results 
are discussed further in Section 17721 



7. DISCUSSION 
7.1. Outflows 

7.1.1. Comparison with Previous Studies 

Since the purpose of this study is to determine whether starbursts or AGN are the pri- 
mary mechanism behind galactic outflows, it is important to compare the results determined 
here to those of RVS05c for the ULIRG Seyferts - galaxies with both starburst and AGN. 
As seen in Table HI the detection rates for the IR-faint Seyfert 2s and Seyfert Is are only 
18% ± 9% and 6% ± 6%, respectively, but rates were as high as 45% ± 11% and 50% ± 
20% for the IR- luminous Seyfert 2s and Is (RVS05c). Our detection rates are even lower in 
comparison to non-Seyfert ULIRGs, since rates of 80% ± 7% and 46% ± 13% were measured 
for low-z and high-z ULIRGs in RVS05b, and a rate of 83% was reported in Martin (2005). 
Poststarburst galaxies have also b een found to have high outflow detection rates on the 



order of 70% (ITremonti et al.l 120071 ). The measured outflow detection rates in the IR-faint 
Seyfert galaxies are therefore considerably smaller than those of pure starbursting galaxies, 
poststarbursts, and AGN + starburst composites. 

Previous observations have indicated that the outflow detection rate increases with 



infrared luminosity (RVS05b. ISato et al.l 120091 ) . Figure M shows the outflow detection rates 



of the IR-faint Seyferts in our sample as well as the the IR-luminous Seyferts and starburst 
ULIRGs and LIRGs of RVS05b and RVS05c as a function of Lfir- The increase in outflow 
detection rate with Lfir is apparent, despite the limited range in Lfir- As Lfir is correlated 
with the star formation rate, this increase in the detection rate with Lfir suggests that star 
formation in all of these systems is the main driver of neutral outflows detected in Na I. 
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Geometry could also be playing a role, since AGN winds with smaller opening angles (higher 
collimation) would reduce the detection rate (see Section 17.1.21) . 

Histograms showing the distributions of all (negative and positive) velocities in our 
current data as well as those of RVS05b and RVS05c provide another point of comparison 
(Figure [9]). The left panel compares the velocities of IR-faint Seyfert galaxy components 
from this study to the Seyfert and starburst ULIRGs & LIRGs of RVS05b and RVS05c. We 
can safely rule out rotation as being a primary cause for the motion of the Na I D gas based 
on the fact that we do not observe a symmetry about zero velocity in these histograms. Were 
rotation the dominant gas motion, we should observe roughly equal amounts of blueshifted 
and redshifted gas in each object and in the overall distributions shown in Figure 8, assuming 
the Na I D gas is distributed more or less symmetrically around the center of rotation and 
is unaffected by severe differential obscuration. Negative velocities are found in a much 
higher percentage of IR-luminous objects (Seyferts of RVS05c and starbursts of RVS05b) 
than IR-faint objects, again indicating that outflows are both stronger and more frequent 
on average in objects with high SFR. The right panel in Figure M combines together the 
IR-faint and IR-luminous Seyferts from this study and RVS05c and compares them to the 
starbursts of RVS05b. In terms of the negative velocities, the results for Seyfert 1 and Seyfert 
2 galaxies are quite different, with Seyfert 2 galaxies showing similar outflow percentages to 
those of the ULIRGs & LIRGs. This suggests a physical connection between the mechanisms 
that drive outflows in Seyfert 2s and starbursting galaxies and a physical difference between 
the mechanisms that cause outflows in Seyfert Is and Seyfert 2s. Interestingly, studies of 
ionized gas out flows have not found such a dichotomy in velocities between Seyfert Is and 



Seyfert 2s (e.g., iRuiz et al.ll200ll ; IVeilleux et al.ll2005l : iRuiz et al.ll2005l : iDas et al.ll2007l . and 



references therein). This suggests that the neutral gas probed by our Na I D observations is 
not kinematically related to the ionized material of the NLR. This difference is even more 
obvious when we also consider our results on positive (inflow) velocities (Section I7.2.ip . 

Rigorous statistical tests confirm these results. Kolmogorov-Smirnov (K-S) and Kuiper 
tests were both used, as the K-S test has an inherent bias in terms of differences in the mean 
and the Kuiper test does not. Low values reported by both of these tests can help rule out two 
sets of data having the same parent distribution. Results of these tests are listed in Table [3 
The small values of P(null) for Av and Av max when comparing IR-faint and IR-luminous 
Seyferts suggest that they do not come from the same parent distribution. A more significant 
comparison comes from combining together the results for the IR-faint Seyfert galaxies from 
this study and those for the IR-luminous Seyfert galaxies from RVS05c, and comparing 
them with those for the starburst ULIRGs and LIRGs from RVS05b. The same statistical 
tests were performed on these distributions, first using all velocities, then the outflowing 
components only. The results are listed in Tables [8] and [9], respectively. The results in Table 
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[8] indicate that Seyferts and starbursts do not share the same parent velocity distributions; 
we return to this point in Section 17.21 When considering only the outflowing components, 
all comparisons show low probability of originating from the same parent distribution except 
when the Seyfert 2 galaxies are compared with the starburst ULIRGs & LIRGs. In that case, 
both the K-S and Kuiper tests return large P(null) values. This confirms quantitatively that 
the outflows in these two classes of objects may arise from the same physical process. 

Figure [TO] displays plots showing Av max as a function of far-IR luminosity (correlated 
with star formation rate) and of galactic circular velocity (correlated with galactic mass) for 
Seyferts from this study, IR-luminous Seyferts from RVS05c, starburst LIRGs & ULIRGs 
from RVS05b, and four starburst dwarf galaxies from Schwartz & Martin (2004). The four 
dwarfs were added in order to see if low-mass galaxies follow the same trends as high-mass 
systems. The IR-faint Seyfert 2 galaxies seem to follow the same trends as the IR-bright 
galaxies, with /S.v max increasing with both SFR and galactic mass, while the outflow velocities 
of all Seyfert 1 galaxies lie above these trends. This again suggests a fundamental difference 
in the way the winds in IR-faint/bright Seyfert 1 galaxies are powered compared with the 
winds in IR-faint/bright Seyfert 2 and starburst galaxies. We return to this issue in Section 

En 



7.1.2. Outflow Dynamics 

Another useful comparison between IR-faint and IR-luminous Seyferts is to look at 
the dynamical properties in addition to the kinematics. Using calculated covering fraction, 
column density, and velocities, we can estimate the mass, momentum, and kinetic energy 
of the neutral gas phase of the ISM being probed by those winds. We follow the method 
outlined in the original study (RVS05b), which made the assumption that the outflows are 
spherically symmetric mass-conserving free winds, with a velocity and instantaneous mass 
outflow rate which do not depend on radius within the wind and which are zero outside 
the wind. This method also assumed that the wind is a thin shell with a uniform radius 
of 5 kpc. This value was based on actual spatial measurements in some of the objects of 
RVS05c, but such spatial measurements are not available for our IR-faint Seyfert galaxies. 
In the present study, we use a 5 kpc radius for both the Seyfert 2 and Seyfert 1 galaxies. 
This is different from the radius of 10 pc that was used for the Seyfert 1 galaxies in the 
original study (RVS05c); the radius in IR-faint Seyfert Is was chosen to be the same as for 
IR-faint Seyfert 2s to facilitate comparisons within the IR-faint sample, and to that end, 
the IR-luminous Seyfert 1 dynamical values from RVS05c have been scaled up to a 5 kpc 
radius as well. These results should be considered order-of-magnitude estimates since they 
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are based on a number of largely unproven assumptions. We have not listed all results here; 
selected results can be seen in Figure [TTJ 

For the objects in the present study, we assume a modest value of 0.3 for the large 
scale opening angle (Cq), which is the typical value for local disk winds (jVeilleux et al. 



20051 ). Following the method of RVS05b, the covering fraction is used to parameterize the 
dumpiness of the wind, or it may reflect the global solid angle subtended by the wind when 
viewed from the galactic center. Using the value of (Cf) as listed in Table IU this yields 
a global covering factor of Q ~ 0.1 for both Seyfert Is and Seyfert 2s. These values are 
not well constrained, especially for the Seyfert Is since only one outflow was detected in 
this type of Seyfert. The Q value is inconsistent with values of Q ~ 0.5-1.0 calculated by 
Crenshaw et al. (2003b) from UV absorption lines in local Seyfert Is. This low Q value could 
imply that the winds we are seeing are collimated rather than wide-angle, and thus our 
lack of outflow detection in some of these objects may be due to their orientation relative 
to us rather than a complete lack of neutral gas outflow. The influence of host galaxy 
inclination on column density was explored using extinction-corrected axis ratios from the 
de Vaucouleurs Third Reference Catalog (?), but no conclusion could be drawn due to 
small-number statistics. Additionally, these Q values will be low in comparison to values 
from other studies; there is contribution here from the background galaxy, and many of 
the outflows we have found are small-scale. Comparisons with the results of UV absorption 
line studies should be considered with caution, though, since UV absorbers are typically 
of much higher ionization (N V, C IV) and likely located much closer to the AGN than 
Na I D absorbing material. The Na I D absorbers may also be affected more strongly by 
host galaxy contamination than the UV absorbers. Blueshifted absorption lines have also 
been detected in X-ray spectra for a number of Seyfert galaxies, including Mrk 6, with column 
densities of si gnificantly higher order of magnitude than we have found here f or Na I D in the 
optical (e.g., 



Mathur et al. 



19951 : iFeldmeier et all Il999l ; Monier et all l200ll ; iKroneold et al. 



20031 ; IVaughan fc Uttleyll2008l ). Again, we caution comparing these results with our own, as 
it has been noted that these X-ray lines are much higher ionization (O VII, O VIII) and are 
often intrinsically related to the aforementioned UV absorption lines. 

One particularly interesting dynamical quantity is the mass outflow rate, since galactic 
outflows may contribu te to the IGM enrichment and are possibly a quencher of star formation 
( iTremonti et al.l 120071 ). Plots of dM/dt for the IR-faint Seyferts, IR-luminous Seyferts, and 
starbursting LIRGs and ULIRGS from RVS05b and RVS05c, and the dwarfs of Schwartz & 
Martin (all calculations are based on a absorber radius of 5 kpc) are presented in Figure [TTJ 
We see a general trend of mass outflow rate increasing with both Lfir and galactic mass. 
There is a considerable difference in the mean dM/dt rates between IR-luminous Seyfert Is 
and IR-luminous Seyfert 2s. The momentum and energy, as well as outflow rates for those 
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quantities, are also significantly higher for IR-luminous Seyfert Is than IR-luminous Seyfert 
2s and all IR-faint Seyferts. The overall good agreement in detection rate and kinematics 
between IR-faint and IR-bright Seyfert 2s and starburst ULIRGs & LIRGs (Figures ISlfTTl) 
suggests that the outflows in all of these objects are powered by star formation, while the 
marked differences when IR-faint Seyfert 1 galaxies are considered suggest that the AGN 
plays an important role in driving the (generally high velocity) outflows in Seyfert 1 galax- 
ies. Escape fractions were not calculated here since too few objects have outflow velocity 
components. 

We can also look at energy outflow rates in order to determine what role these outflows 
could play in galactic feedback. The rates that we have calculated can be found in Table [31 
though it must be cautioned that these numbers are a function of our uncertain Q value, as 
well as our assumed absorber radius of 5 kpc. For the Seyfert 2 outflows, the average energy 
outflow rate was found to be ~ 10 41 - 1 (g^) ergs s" 1 , or ~ 10 7 (^j) (g^) L & . In 
comparison to the average bolometric luminosity for these objects (~ 10 10 L & , taken from 



Woo fc Urryl |2002| ) . the energy outflow rates are only ~1% of the host galaxy luminosity. 
This indicates that outflow energetics in our Seyfert 2 galaxies are not strong enough to 
play a large role in galactic feedback. The findings of Schlesinger et al. (2009) for Mrk 573 
are in agreement with our conclusion. For the only Seyfert 1 galaxy that we have found 
to have an outflow, Mrk 6, the energy outflow rate of its lower velocity component (10 4L6 
(<fr) ( 5 k c ) er § s s_1 ) ^ s a S a i n ~1% °f the host bolometric luminosity. This indicates that 
this Seyfert 1 outflow is also not energetic enough to play a vital role in galact i c feedback. 



and t hus previous findings are consistent with our results (IKrongold et al.l 120071 ; IStoll et al. 



20091 ). If we look at the higher velocity outflow found in Mrk 6, its energy outflow rate 
(10 42 ' 9 (qy) ( 5 ^ - ) ergs s _1 ) is ~5% of the host bolometric luminosity, which is higher than 
for the Seyfert 2 galaxies but again likely not strong enough to influence the evolution of its 
environment. However, one should note that these values are all highly dependent on the 
calculated value of Q and assumed value of r and are thus uncertain in comparison to our 
measured velocities. 



7.2. Inflows 

7.2.1. Comparison with Previous Studies 

An unexpected result of the present study is the high detection rate of inflows in IR-faint 
Seyferts (39% ± 12% for Seyfert Is, 35% ± 11% for Seyfert 2s). In contrast, only -15% 
of the IR-luminous objects in RVS05b and RVS05c showed redshifted Na I D absorption. 
This difference is clearly seen in the left panel of Figure [91 Interestingly, a recent search for 
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outflows in the AEGIS da tabase has also revealed an excess of inflows among AGN-powered 
systems ( ISato et al.ll2009l ). Inflow has been observed in at least one Seyfert 1 galaxy, NGC 
5548, using ionized gas detected in the UV, though detections have not been reported for such 
a large number of objects as we have found here (IMathur et al.lll999l ). There has been one 
tentative observation of redshifted X-ray absorption in a Seyfert 1 galaxy, but the authors 
caution that the sig nificance of the absorption line they have measured is highly uncertain 
rtDadina et al1l2005h . 



When we examine the inflow data on the IR-faint galaxies alone, we find no significant 
correlation between Av max and galactic mass (left panel in Figure [12]). Neither is there 
any obvious trend with Seyfert type or far- infrared luminosity (right panel of Figure [T2l . in 
contrast to the trends seen for the outflowing gas. Additionally, we see inflows in nearly the 
same fraction of Seyfert Is and Seyfert 2s, indicating no particular trend with Seyfert type. 
As for the detected outflows, the influence of host galaxy inclination on column density was 
explored using extinction-corrected axis ratios (?), but again, no conclusion could be drawn 
due to small-number statistics. 

Next, the velocity distributions of the inflowing components for the IR-faint and IR- 
luminous Seyfert 1 galaxies were combined together and compared with the combined distri- 
bution of IR-faint and IR-luminous Seyfert 2 galaxies. K-S and Kuiper tests were performed 
on Av, Av max , and Doppler parameter, following the same procedure as in Section I7TTI and 
the results are listed in Table (TUl They confirm the lack of obvious differences in the inflow 
properties between the two types of Seyfert galaxies. 



7.2.2. Search for Connection between Inflows and Nuclear Structures 



The Na I absorption infall velocities often extend to relatively high values so we fa- 
vor a nuclear location for this gas rather than a galactic origin (e.g. galactic fountains, 
Ferrara fc Ricottill2006l ). Nuclear accretion of cool gas like Na I with T ~ 100 K ca n provide 



fuel not only for star formation but also for nuclear activity (jStruve et al.ll2008l ). Various 
mechanisms have been proposed to help reduce the large angular momen tum of the gas in 
the nuclei of galaxies. The se include nuclear (< 1 kpc) bars and spirals (IPogge fc Martini 



2002 



Martini et all l2003al fbL ?) and gravitational interactions with neighboring galaxies 



( iCrenshaw et al.ll2003al ). Nuclear bars are thought to form when a large galactic bar forces 
gas inwards, creating a gaseous disk, and instability cau ses formation of a small gas bar near 
the nucleus ( IShlosman et all Il989l ; iMartini et al.l 1200 ll ) . Stellar bars could also accomplish 
the same thing. B oth bar types are cap able of removing angular momentum from gas rotating 
near the nucleus (iMartini et all 1200 ll ). Nuclear spirals have been proposed as another AGN 
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fueling mechanism since shock fronts that occur at their boundaries can take away angular 



mome ntum from local gas and cause material to fall in towards the black hole (IMartini et al. 



2003bl ). However, these nuclear structures do not necessarily l ead to AGN fueling sinc e they 
are present in a equally large fraction of non-active galaxies ( IPogge &: Martini! |2002| ). Star 
formation may occur in these objects and disrupt AGN- fueling inflows (?). 

We revisit this issue here by looking for the presence of nuclear spirals or bars in the 
IR-faint Seyfert galaxies with inflow signature. We have compared our results to nuclear 
structure survey s done by three different groups (IMalkan et al.1 Il998l ; IMartini et al.1 l2003al ; 
Deo et al.ll2006l ). Nuclear dust structures in objects in common with these studies are classi- 
fied into five distinct morphological categories (see Column (11) of Tabled]): irregular dust, 
dust filaments, nuclear dust ring, nuclear dust bar, and nuclear dust spiral. Of the thirteen 
objects which show inflows, five show evidence for nuclear dust spirals, bars, or rings, one 
shows evidence for dust filaments, one for irregular dust, and six show no sign of nuclear 
dust structure. The significant fraction of objects that show both inflows and nuclear dust 
spirals, bars, rings, and/or filaments lends credence to the idea of a connection between 
morphology and kinematics. Of the seven objects with irregular dust or no dust, five have 
nearby companion galaxies, so tidal forces due to interactions with these companions ma y 
cause AGN fueling (IRafanelli et al.lll995t iHenriksen Cousineaulll999l ; ISmith et al.1 120071 ). 
This leaves only two objects with inflow which do not show nuclear structure or a compan- 
ion: Akn 202 and Mrk 1018. There are also two objects in our study with nuclear dust 
structure that show measurable outflow rather than inflow (NGC 3786 and NGC 7319), and 
thus whether we measure outflow rather than inflow may be a consequence of our line of 
sight to the nucleus (inflow and outflow may be occurring in different planes, or in the same 
plane but over a different range of azimuthal angles), rather than a lack of inflow in the 
object. Measurements of the line-of-sight velocity field with a resolution of ~10s of pc will 
be needed to disentangle the geometry of the inflow s /outflows detected in our data (e.g., the 
study of NGC 1097 by jvan de Ven fc Fathill2009[ ). 



7.2.3. Inflow Dynamics 



The same method used to calculate mass, momentum, and kinetic energy for outflows 
(Section 17.1. 2p was used for inflows, but the characteristic absorber radius was reduced to 1 
kpc, a rough upper limit to the scale of the circumnuclear str uctures (nuclear bars /spirals) 
believed to be responsible for feeding the AGN (Section |7.2.2[ IMartini et al.ll2003a ). Again, 
these dynamical quantities are rather uncertain, but are calculated to find out at least 
roughly how these inflows compare with the mass accretion rates necessary to power the 
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AGN in these systems. Table |5] lists the mass, momentum, and kinetic energy calculated for 
all objects with measured inflows. The mass accretion rate for these inflows ranges from just 



under 1 



1 kpc< 



M yr 1 to just under 5 



M yr" 1 . For comparison, the mass accretion 
Lhn]/c 2 r], where L hoi is the bolometric luminosity 



rate needed to power an AGN is M = ^boi/^ 'h Wli ere lJboi 
of the AGN, M is the mass accretion rate, c is the speed of light, and rj is an efficiency 
factor dictating how much of the rest mass of the material being accreted is turned into 
radiation. If we take rj to be ~ 0.1 ( Riffel et al. 2008h and allow the bolometric luminosity to 



be 



10 ergs s , typical for Seyfert galaxies (jPadovani fc Rafanellilll988l ; iCrenshaw et al. 



2003bl ). then we find that a mass accretion rate of M ~ 10 2 M Q yr 1 is required. Even 



with our rough order of magnitude estimates, the mass accretion rates of all of our observed 
Na I D inflows are well above the amount necessary to power the AGN. Thus the inflows 
that we are measuring carry enough material to fuel the AGN in these objects, even if only 
~1% of this material makes its way down to the AGN. The total infalling mass of ~ 10 7 
M Q , estimated fro m our data, is enough to su stain nuclear activity over typical AGN 



1 kpc ' 



lifetimes (~ 10' 



10* 



yrs; 



Mouri fc Taniguchill2002l : ?). 



8. SUMMARY 

The main results from our study of Na I D absorption in infrared-faint Seyfert 1 and 
Seyfert 2 galaxies can be summarized as follows: 

• Outflow Detection Rates and Kinematics: The rates of detection of outflows in IR-faint 
Seyfert 1 and 2 galaxies are lower than previously found in IR-bright Seyferts. Outflows 
were found in ~18% of IR-faint Seyfert 2s in our sample, compared with ~45% among 
the Seyfert 2 ULIRGs of RVS05c. Only one out of 18 Seyfert 1 galaxies in our sample 
shows evidence for a wind, far lower than the ~50% reported for Seyfert 1 ULIRGs 
in RVS05c. Interestingly, the outflow detection rate and velocities of IR-faint Seyferts 
follow the same trends with infrared luminosity and galaxy mass as IR-bright systems. 
The outflow kinematics of Seyfert 2 galaxies resemble those of starburst galaxies, while 
the outflow velocities in Seyfert 1 galaxies are significantly larger. These results suggest 
that the AGN is contributing to the neutral outflows in Seyfert 1 systems, while the 
starburst is the main driver of the outflows in all Seyfert 2 galaxies. Differences in 
wind angular extent (e.g. AGN-driven outflows in Seyfert Is are more collimated than 
starburst- driven outflows) may also explain some of these results. 

• Outflow Dynamical Estimates: The mass, momentum, and kinetic energy of the mate- 
rial involved in these outflows were estimated assuming a constant characteristic radius 
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of 5 kpc. The dynamical properties of the outflows in Seyfert 2 galaxies are similar 
to those of the starburst ULIRGs & LIRGs, but differ significantly from those of the 
Seyfert Is. This again suggests a fundamental physical difference between the outflows 
of Seyfert Is and those in the other objects. An attempt to determine the influence of 
host galaxy inclination on these outflows was inconclusive. Measured energetic rates 
do not appear large enough to play a significant role in galactic feedback, but these 
values are uncertain. 

• Inflow Detection Rates and Kinematics: A striking result of this study is the high rate 
of detection of spatially-unresolved redshifted Na I D absorption, which we interpret 
as nuclear inflows (39% and 35% inflow detection rates for Seyfert 2s and Is), with 
maximum velocities (321 and 291 km s _1 on average). This is evidence for the existence 
of some mechanism capable of removing angular momentum from the circumnuclear 
gas in these objects. Nuclear bars and spirals, as well as interactions with nearby 
galaxies, may play a role in this process. 

• Inflow Dynamical Estimates: Mass, momentum, and kinetic energy were estimated for 
the inflows, using a characteristic radius of 1 kpc, consistent with the observed sizes 
of nuclear bars and spirals in these systems. While these estimates are uncertain, we 
find that the total infalling mass and infalling mass rates are more than enough to 
power the AGN in these systems for typical AGN lifetimes. As with the outflows, an 
attempt was made to determine the influence of inclination on column density, but no 
conclusion could be drawn. 
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Table 1. Galaxy properties and observing logs 



Name 


z 




£fir 


SFR 


V c 


w eq 


Run 




PA 


Structure 


Refs 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Seyfert 2s 


Akn 79 


0.01743 








231 


3.43 


2005 Sep 05 


4500 







1 


Arp 107A 


0.03463 


10.63 


10.09 


2.14 


308 


2.20 


2006 Apr 02 


5400 







1 


Mrk 348 


0.01516 


10.40 


9.89 


1.33 


191 


1.53 


2006 Nov 23 


4800 







1 


Mrk 622 


0.02347 


10.72 


10.26 


3.15 


301 


2.03 


2006 Apr 02 


5400 





I 


1,4,7 


Mrk 686 


0.01420 


10.04 


9.64 


0.76 


558 


2.98 


2006 Mar 31 


4500 


90 


S 


1,5,6 


Mrk 1157 


0.01510 


10.41 


10.13 


2.32 


585 


1.99 


2005 Sep 04 


3540 





R 


1,4,7 


NGC 1358 


0.01351 


10.35 


9.37 


0.41 


299 


3.73 


2006 Nov 21 


4800 







1 


NGC 1667 


0.01517 a 


10.97 


10.66 


7.84 


580 


3.99 


2006 Nov 21 


4800 





S 


1,4,6 


NGC 3362 


0.02767 








358 


1.92 


2006 Mar 31 


4500 


90 


F 


1,7 


NGC 3786 


0.00903 








367 


3.11 


2006 Apr 02 


5400 


22 


R,S 


1,6,8 


NGC 4388 


0.00839 


10.41 


10.34 


3.78 


414 


2.27 


2006 Apr 02 


5400 





I 


1,5,6 


NGC 5252 


0.02308 


10.39 


9.84 


1.18 


340 


3.72 


2006 Apr 02 


5400 





R,S 


1,6,7 


NGC 5728 


0.00932 a 


10.60 


10.32 


3.63 


434 


2.25 


2006 Apr 02 


3600 







2,4 


NGC 7319 


0.02236 


10.56 


10.21 


2.83 


210 


6.61 


2005 Sep 08 


3240 





F 


1,7 


NGC 7672 


0.01348 


9.91 


9.64 


0.76 


363 


1.39 


2006 Nov 21 


7200 







1,4 


NGC 7682 


0.01707 








412 


2.06 


2006 Nov 21 


6000 





S 


1,6 


UGC 3995 


0.01575 


10.32 


9.76 


0.99 


250 


3.25 


2006 Mar 31 


4800 


90 




1,4 


Seyfert Is 


Akn 202 


0.02872 b 


10.71 


10.15 


2.43 


610 


0.44 


2006 Nov 23 


6000 







4 


Akn 564 


0.02468 a 


10.75 


10.14 


2.40 


443 


0.00 


2005 Sep 08 


2040 







2,5 


Mrk 6 


0.01951 a 


10.59 


10.09 


2.14 


593 


1.08 


2006 Nov 21 


6000 





I 


1,4,7 


Mrk 10 


0.02925 


10.78 


10.38 


4.14 


620 


2.18 


2006 Nov 21 


3600 







1 


Mrk 79 


0.02221 a 


10.85 


10.33 


3.68 


356 


0.39 


2006 Nov 22 


4800 





S 


1,4,8 


Mrk 110 


0.03513 b 








147 


0.00 


2006 Nov 22 


3600 









Mrk 352 


0.01486 a 








300 


0.60 


2006 Nov 22 


7200 







2 


Mrk 359 


0.01694 


10.35 


9.96 


1.59 


235 


0.55 


2005 Sep 05 


2700 





s 


2,5,8 


Mrk 382 


0.03348 b 


10.72 


10.00 


1.73 




0.00 


2006 Nov 23 


4800 





B,S 


3,7,8 


Mrk 477 


0.03744 b 


11.14 


10.72 


9.01 




0.00 


2006 Mar 31 


3926 


90 




4 


Mrk 506 


0.04303 a 


10.69 


10.12 


2.30 


170 


2.96 


2005 Sep 08 


4830 


27 




1 


Mrk 595 


0.02739 b 


10.64 


10.28 


3.31 


444 


1.41 


2006 Nov 23 


6000 





s 


3,8 


Mrk 1018 


0.04263 








414 


1.58 


2006 Nov 23 


6000 







1 


Mrk 1126 


0.01057 








540 


1.70 


2006 Nov 22 


7200 





s 


1,8 


NGC 788 


0.01350 


10.04 


9.38 


0.41 


226 


3.20 


2005 Sep 04 


2700 





s 


1,6 


NGC 1019 


0.02460 a 


10.43 


9.97 


1.62 


178 


2.63 


2006 Nov 22 


6000 





B,R,S 


2,7,8 


NGC 7603 


0.02956 


10.78 


10.41 


4.43 


313 


0.45 


2006 Nov 23 


7200 





R,S 


1,4,8 



Table 1 — Continued 



Name 


z 




£fir 


SFR 


v c 


w eq 


Run 


^exp 


PA 


Structure 


Refs 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


UGC 3223 


0.01567 


10.37 


10.09 


2.13 


263 


1.66 


2006 Nov 22 


6000 





S 


1,4,8 



Note. — Col. (2): Heliocentric redshift. All redshifts are based on stellar measurements except (a) HI 21-cm 
measurements, or (b) measured from the emission lines in our data. Col. (3): Infrared luminosity, in logarithmic 
units of L®. Col. (4): Far- infrared luminosity, in logarithmic units of L Q (see Section [2]2]). Col. (5): Star formation 
rate, computed from the far-infrared luminosity (Section I2.2|) . Col. (6): Circular velocity, equal to V2a, v ro t, or 
the quadratic combination if both are available. Col. (7): Rest-frame equivalent width of Na ID as computed 
from our model fits. Col. (8): Observing dates (Section [3]). Instrument used was the R-C Spectrograph on the 
KPNO 4m. Col. (9): Total exposure time in seconds. Col. (10): Slit position angle. Col. (11): Letters indicate 
dust structure around nucleus (see Section !?. 2. 2| : (B) nuclear dust bar, (F) dust filaments, (I) irregular dust, (R) 
nuclear dust ring, (S) nuclear dust spiral. Col. (12): Reference. 



References. — m lNelson fc Whittlelll995l: (2) ISpringob et alJliool (3) IWhittlelll992at (4) IRAS Faint Sou rce 
Catalog; (5) IRAS Point Source Catalog; (6) iMartini et al.ll2003al: (7) iMalkan et al.lll998k (8) beo et al.ll2006[ 
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Table 2. Properties of individual velocity components 



Ai, c Aw b n, c C f TV(Nal) JV(H) 

Name (A) (km s _1 ) (km s _1 ) (cm ~ 2 ) (cm ~ 2 ) 

(1) (2) (3) (4) (5) (6) ' (7) ' (8) 



Seyfert 2s 



Akn 79 


5999 


.47 


-44 ± 5 


85 ± 


19 


1. 


.07 


+0.50 
— 0.21 





.12 


+() 

— 0. 


06 
.02 


13 


.44 


+0. 
— 0. 


28 
.20 


21 


.08 


+0.43 
— 0.43 






5999 


.86 


-25 ± 7 


280 ± 


22 


0. 


.27 


+0.04 
— 04 





.48 


+0 
— 


07 

06 


13 


.42 


+0. 
— 


12 
12 


21 


.06 


+0.18 

— 18 


Arp 107A 


6102 


.90 


54 ± 24 


323 ± 


59 


0. 


.07 


+0'05 

Q Ql 





.92 


+0^ 

o 


04 
01 


12. 


.94 


+0. 
() 


25 
9 


20 


.58 


+0'25 
q og 


Mrk 348 


5988 


.26 


65 ± 4 


261 ± 


25 


0. 


.06 


+0.01 
-0.01 


1 


.00 


+0. 
-0 


01 
.01 


12. 


.74 


+0 
-0 


08 
06 


20 


.38 


+0.08 
-0.06 


Mrk 622 


6033 


.19 


-138 ± 10 


416 ± 


62 


0. 


.09 


+0.04 
-0.01 





.55 


+0. 
-0 


23 
.02 


13 


.13 


+0 
-0 


16 

.07 


20 


.77 


+0.16 
-0.07 


Mrk 686 


5981. 


.63 


16 ± 5 


287 ± 


34 


0. 


.10 


+0.04 
-0.01 


1 


.00 


+0 
-0 


01 
.01 


13 


.05 


+0. 
-0. 


25 
11 


20 


.69 


+0.31 
-0.17 


Mrk 1157 


5986 


.67 


3 ± 2 


211 ± 


20 


0. 


.83 


+0.17 
-0.13 





.19 


+0 
-0 


04 
.03 


13 


.79 


+0. 
-0 


17 

13 


21 


.43 


+0.23 
-0.20 


NGC 


1358 


5977. 


.92 


35 ± 6 


320 ± 


19 


0. 


.43 


+0.05 
-0.05 





.37 


+0 
-0 


04 
.04 


13. 


.69 


+0 
-0 


08 
.08 


21 


.33 


+0.15 
-0.15 


NGC 


1667 


5986 


.79 


-12 ± 84 


343 ± 


19 


0. 


.32 


+0.03 
-0.03 





.46 


+0 
-0 


05 
.05 


13 


.59 


+0. 
-0 


06 

.06 


21 


.23 


+0.13 
-0.13 


NGC 


3362 


6062 


.43 


84 ± 38 


323 ± 


92 


0. 


.05 


+0.07 
-0.01 


1 


.00 


+0 
-0 


01 
.01 


12. 


.83 


+0. 
-0 


39 
16 


20 


.47 


+0.39 
-0.16 


NGC 


3786 


5947. 


.88 


-148 ± 5 


34 ± 3 


5. 


.00 


+0.36 
-2.00 





.10 


+0 
-0 


01 

.04 


>13 


.76 


+0. 
-0 


05 
.23 


>21 


.40 


+0.05 
-0.23 






5951 


.19 


19 ± 2 


239 ± 14 


0. 


.34 


+0.04 
-0.03 





.44 


+0 
-0 


05 
.04 


>13 


,45 


+0 
-0. 


03 
21 


>21 


.09 


+0.03 
-0.21 


NGC 


4388 


5948 


.05 


51 ± 17 


125 ± 5 


0. 


.98 


+0.07 
-0.07 





.32 


+0 
-0 


02 
.02 


13. 


.60 


+0. 
-0 


04 
.04 


21 


.24 


+0.04 
-0.04 


NGC 


5252 


6034 


.47 


39 ± 10 


375 ± 23 


0. 


.10 


+0.02 
-0.01 


1 


.00 


+0 
-0 


01 
.01 


13 


.14 


+0. 
-0 


15 
.07 


20 


.78 


+0.21 
-0.13 


NGC 


5728 


5952 


.30 


-11 ± 1 


262 ± 26 


0. 


.27 


+0.06 
-0.03 





.37 


+0. 
-0 


08 
.04 


13 


.40 


+0. 
-0 


17 
.13 


21 


.04 


+0.23 
-0.19 


NGC 


7319 


6026 


.74 


-133 ± 19 


100 ± 1 


3. 


.00 


+0.24 
-2.12 





.20 


+0 
-0 


02 
14 


13. 


34 


+0. 
-0. 


04 
.10 


20 


.98 


+0.04 
-0.10 






6026 


.82 


-130 ± 1 


450 ± 1 


0. 


.46 


+0.06 
-0.06 





.49 


+0. 

-0 


06 

.06 


13. 


.81 


+0 
-0 


05 
12 


21 


.45 


+0.05 
-0.12 


NGC 


7672 


5976 


.55 


-25 ± 24 


98 ± 6 


1. 


.12 


+0.17 
-0.07 





.24 


+0 
-0 


04 
.01 


13 


.53 


+0 
-0 


11 

.06 


21 


.17 


+0.18 
-0.03 


NGC 


7682 


6000 


.21 


99 ± 3 


228 ± 21 


0. 


.49 


+0.09 
-0.08 





.25 


+0. 
-0. 


05 
.04 


13. 


.59 


+0. 
-0 


08 
.09 


21 


.23 


+0.08 
-0.09 


UGC 


3995 


5992 


.99 


127 ± 2 


268 ± 15 


0. 


.13 


+0.02 
-0.01 


1 


.00 


+0 
-0 


01 
.01 


13 


.09 


+0. 
-0 


07 

03 


20 


.73 


+0.07 
-0.03 


Seyfert Is 


Akn : 


202 


6069 


.22 


113 ± 11 


129 ± 32 


1. 


.02 


+0.23 
-0.05 





.05 


+0 
-0 


01 
.01 


13 


.67 


+0 
-0 


.12 
.13 


21 


.31 


+0.12 
-0.13 


Mrk i 


5 


6008 


.03 


-229 ± 1 


25 ± 1 


1. 


.46 


+0.03 
-0.03 





.48 


+0 
-0 


01 
.01 


13 


.12 


+0 
-0 


17 

29 


20 


.76 


+0.17 
-0.29 






5992 


.12 


-1024 ± 1 


16 ± 


11 


4. 


.75 


+0.23 
-0.89 





.06 


+0. 
-0 


01 
.01 


13 


.44 


+0 
-0 


.20 
.31 


21 


.08 


+0.20 
-0.31 


Mrk 


10 


6070 


.59 


26 ± 7 


305 ± 


34 


0. 


.17 


+0.06 
-0.02 





.45 


+0. 
-0 


15 
.04 


13 


.28 


+0. 
-0 


23 
11 


20 


.92 


+0.30 
-0.17 


Mrk 


79 


6028 


.25 


-15 ± 7 


112 ± 


23 


0. 


.92 


+0.32 
-0.31 





.06 


+0 
-0 


06 
.02 


13 


.57 


+0 
-0 


25 
.24 


21 


.21 


+0.31 
-0.30 


Mrk : 


352 


5986 


.09 


45 ± 8 


111 ± 


25 


0. 


.16 


+0.04 
-0.01 





.16 


+0 
-0 


03 
.01 


12. 


.82 


+0 
-0 


07 
.42 


20 


.46 


+0.07 
-0.42 


Mrk : 


359 


5997 


.12 


-17 ± 11 


147 ± 


28 


0. 


.21 


+0.12 
-0.03 





.20 


+0 
-0 


12 
.03 


13 


.03 


+0. 
-0 


30 
.17 


20 


.67 


+0.36 
-0.23 


Mrk - 


506 


6154 


.09 


135 ± 7 


296 ± 


32 


0. 


.17 


+0.05 
-0.01 





.64 


+0 
-0 


18 

.05 


13 


.26 


+0 
-0 


12 

.06 


20 


.90 


+0.12 
-0.06 


Mrk . 


595 


6059 


.26 


8 ± 5 


230 ± 


21 


0. 


.15 


+0.03 
-0.01 





.44 


+0. 
-0 


08 
03 


13 


.09 


+0. 
-0 


15 
05 


20 


.73 


+0.22 
-0.12 


Mrk 


1018 


6150 


.34 


67 ± 5 


379 ± 


31 


0. 


.12 


+0.03 
-0.01 





.36 


+0. 
-0 


09 
.03 


13 


.22 


+0 
-0 


10 
.05 


20 


.86 


+0.10 
-0.05 


Mrk 


1126 


5962 


.64 


138 ± 3 


204 ± 


11 


0. 


.08 


+0.01 
-0.01 


1 


.00 


+0 
-0 


01 
.01 


12. 


.79 


+0. 
-0 


06 
.03 


20 


.43 


+0.06 
-0.03 


NGC 788 


5977 


.35 


± 3 


253 ± 


61 


0. 


.72 


+0.06 
-0.08 





.28 


+0 
-0 


02 
.03 


13 


.81 


+0. 
-0 


19 
19 


21 


.45 


+0.25 
-0.25 


NGC 1019 


6044 


.02 


68 ± 4 


253 ± 


23 


0. 


.31 


+0.05 
-0.04 





.41 


+0 
-0 


07 

.05 


13. 


.44 


+0 
-0 


07 
.07 


21 


.08 


+0.07 
-0.07 
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Table 2 — Continued 



Name 
(1) 


Ai, c 
(A) 

(2) 


Av 

(km s -1 ) 
(3) 


b 

(km s- 1 ) 
(4) 


n, c 

(5) 


c f 

(6) 


iV(Na I) 
(cm - 2 ) 
(7) 


JV(H) 
(cm -2) 

(8) 


NGC 7603 
UGC 3223 


6073.78 
5989.19 
5991.73 


93 ± 11 
-39 ± 1 
88 ± 1 


497 ± 45 
81 ± 10 
98 ± 22 


o.o6 ±g:8S 
5.00 i?:S 
0.92 t° il 


o.i4 ±8:iS 
o.o8 ±g:S 
o.i6 ±8:S 


13.06 ±0;22 

14.16 t° .ll 
13.50 +°;« 


20.70 +°-^ 
21.80 ±gii? 
21.14 



Note. - Col. (2): Redshiftcd heliocentric wavelength, in vacuum, of the Na I Di A5896 line. 
Col. (3): Velocity relative to systemic. Negative velocities are blueshifted, positive are redshiftcd. 
Components with Av < —50 km s _1 and \Av\ > 26 (Av) are assumed to be outflowing; those with 
Av > 50 km s _1 and |Au| > 26(Av) are assumed to be inflowing . Col. (4): Doppler parameter. Col. 
(5): Central optical depth of the Na I Di A5896 line; the optical depth of the D2 line is twice this 
value. Col. (6): Covering fraction of the gas. Col. (7-8): Logarithm of column density of Na I and H, 
respectively. 



Table 3. Outflow: Individual objects 







M 


dM/dt 


P 


dp/dt 


E 


dE/dt 


Name 


(km S -1 ) 


(M ) 


(M Q yr" 1 ) (dyn s) 


(dyn) 


(ergs) 


(ergs s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Seyfert 2s 


Mrk 622 


-484 


8.76 


1.22 


49.2 


34.2 


56.1 


41.0 


NGC 3786 


-176 


8.64 


1.12 


49.1 


34.1 


56.0 


41.0 


NGC 7319 


-217 


8.52 


0.95 


48.9 


33.9 


55.8 


40.7 




-504 


9.39 


1.85 


49.8 


34.8 


56.7 


41.7 


Seyfert Is 


Mrk 6 


-250 


8.69 


1.36 


49.4 


34.5 


56.4 


41.6 




-1037 


8.08 


1.40 


49.4 


35.2 


57.1 


42.9 



Note. — Col.(2): Maximum velocity in the outflow, Av max = Av - FWHM/2. Col 
.(3): Log of total outflowing mass. Col. (4): Log of mass outflow rate. Col. (5): Log of 
total momentum of outflow. Col. (6): Log of momentum outflow rate. Col. (7): Log of total 
kinetic energy of outflow. Col. (8): Log of kinetic energy outflow rate. Note that all values 
are calculated using fi ~ 0.1. 
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Table 4. Outflow: Average properties 



Quantity 
(1) ' 


IR-Faint Scyfert 2s 
(2) 


IR-Faint Seyfert Is 
(3) 


IR-Lum. Scyfert 2s 
(4)' 


IR-Lum. Seyfert Is 
(5)' 


Number of galaxies 
Detection rate (%) 


17 

18 ± 9 


18 
6 ± 6 


20 
45 ± 11 


6 

50 ± 20 


Galaxy Properties 


z 

log(X Fm /£ Q ) 
SFR (Mq yr- 1 ) 
Av (km s" 1 ) 
Av max (km s _1 ) 
log[iV(Na I D )/cm s~ 2 ] 
log[iV(H)/cm s~ 2 ] 


0.018 ± 0.01 
10.02 ± 0.35 

2.39 ± 1.92 
-137 ± 8 

-345 ± 173 
13.51 ± 0.33 
21.15 ± 0.33 


0.025 ± 0.01 
10.14 ± 0.23 
2.99 ± 2.13 
-627 ± 562 
-643 ± 556 
13.28 ± 0.23 
20.92 ± 0.23 


0.148 ± 0.11 
11.86 ± 0.31 
118 ± 109 
-322 ± 388 
-618 ± 422 
13.5 ± 0.7 
20.9 ± 0.7 


0.150 ± 0.09 
12.13 ± 0.1 

164 ± 75 
-4942 ± 2831 
-5210 ± 4306 
14.5 ± 0.8 
21.8 ± 0.8 


Velocity Component Properties 


T 

b (km s -1 ) 
Cf 


0.76 t° A 7 
250 ± 214 
0.34 


1.01 ±g;5 
21 ± 6 

0.27 t° Q :l 


0.27 ±»:5 

232 ± 182 
0.42 t° 5 2 


0.69 ±J;« 
87 ± 130 
0.67 



Note. - - "IR-Lum." refers to the IR-luminous data of RVS05c. For most quantities we list the mean and la 
dispersions, under the assumption of a Gaussian distribution in the log of the quantity. Statistics for all quantities 
except z, Lfir, and SFR are computed only for outflowing velocity components. Note that the entries under IR-Faint 
Seyfert 2s (Is) corresponds to only three (one) objects. 
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Table 5. Inflow: Individual objects 







M 


dM/dt 


P 


dp/dt 


E 


dE/dt 


Name 


(km s -1 ) 


(M e ) 


(M yr- 1 ) (dyn s) 


(dyn) 


(ergs) 


(ergs s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Seyfert 2s 


Arp 107A 


323 


7.39 


0.15 


48.4 


32.7 


53.9 


39.2 


Mrk 348 


282 


7.22 


0.04 


47.3 


32.7 


53.8 


39.2 


NGC 3362 


352 


7.32 


0.27 


47.6 


33.0 


54.2 


39.6 


NGC 4388 


155 


7.59 


0.31 


47.6 


32.8 


54.0 


39.2 


NGC 7682 


289 


7.48 


0.49 


47.8 


33.3 


54.5 


40.0 


UGC 3995 


350 


7.58 


0.69 


48.0 


33.6 


54.8 


40.4 


Seyfert Is 


Akn 202 


220 


6.89 


-0.04 


47.3 


32.8 


54.0 


39.6 


Mrk 506 


381 


7.55 


0.69 


48.0 


33.6 


54.8 


40.5 


Mrk 1018 


382 


7.27 


0.10 


47.4 


32.7 


53.9 


39.3 


Mrk 1126 


308 


7.27 


0.43 


47.7 


33.4 


54.6 


40.2 


NGC 1019 


279 


7.54 


0.39 


47.7 


33.0 


54.2 


39.6 


NGC 7603 


507 


6.70 


-0.32 


47.0 


32.5 


53.6 


39.1 


UGC 3223 


169 


7.19 


0.15 


47.4 


32.9 


54.1 


39.5 



Note. — All inflow dynamics calculations assume an absorber radius of 1 kpc. Col. (2): 
Maximum velocity of the inflow, Av max = Av + FWHM/2. Col. (3): Log of total inflowing 
mass. Col. (4): Log of mass inflow rate. Col. (5): Log of total momentum of inflow. Col. (6): 
Log of momentum inflow rate. Col. (7): Log of total kinetic energy of inflow. Col. (8): Log 
of kinetic energy inflow rate. 



Table 6. Inflow: Average properties 



Quantity 


IR-Faint Seyfert 2s 


IR-Faint Seyfert Is 


(1) ' 


(2) 


(3) 


Number of galaxies 


17 


18 


Detection rate (%) 


35 ± 11 


39 ± 12 


Galaxy Properties 


z 


0.018 ± 0.01 


0.025 ± 0.01 


log(L F iit/£©) 


10.02 ± 0.35 


10.14 ± 0.23 


SFR (M Q yr- 1 ) 


2.39 ± 1.92 


2.99 ± 2.13 


Ai> (km s -1 ) 


80 ± 29 


100 ± 29 


kVmax (km S _1 ) 


291 ± 37 


321 ± 29 


log[iV(NaI D )/cm s~ 


- 2 ] 13.13 ± 0.38 


13.28 ± 0.29 


log[iV(H)/cm s" 2 ] 


20.77 ± 0.38 


20.92 ± 0.29 



Velocity Component Properties 


T 


n o +o.i3 

U - J -0.06 


o.38 ±g:S 


b (km s- 1 ) 


254 ± 74 


265 ± 140 


Cf 


75 + - 09 

u -'° -0.05 


n so + 02 

u - ,3y -0.02 



Note. — For most quantities we list the mean and la dispersions, un- 
der the assumption of a Gaussian distribution in the log of the quantity. 
Statistics for all quantities except z, Lfir, and SFR are computed only 
for inflowing velocity components. 
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Table 7. Statistical comparisons of kinematic parameters: IR-faint and IR-luminous 

Seyfert galaxies 



Samples 


P(null, K-S) 


_P(null, Kuiper) 


Av 


IR-Faint Seyfert Is vs IR-Faint Seyfert 2s 


0.76 


0.69 


IR-Faint Seyfert Is vs IR-Luminous Seyfert Is 


<0.01 


<0.01 


IR-Faint Seyfert 2s vs IR-Luminous Seyfert 2s 


0.04 


0.07 


IR-Faint Seyferts vs IR-Luminous Seyferts 


<0.01 


<0.01 




IR-Faint Seyfert Is vs IR-Faint Seyfert 2s 


0.84 


0.95 


IR-Faint Seyfert Is vs IR-Luminous Seyfert Is 


<0.01 


<0.01 


IR-Faint Seyfert 2s vs IR-Luminous Seyfert 2s 


<0.01 


0.03 


IR-Faint Seyferts vs IR-Luminous Seyferts 


<0.01 


<0.01 


Doppler Parameter 


IR-Faint Seyfert Is vs IR-Faint Seyfert 2s 


0.33 


0.71 


IR-Faint Seyfert Is vs IR-Luminous Seyfert Is 


0.09 


0.06 


IR-Faint Seyfert 2s vs IR-Luminous Seyfert 2s 


0.28 


0.72 


IR-Faint Seyferts vs IR-Luminous Seyferts 


0.23 


0.41 



Note. — P(null) is the probability that the two listed distributions are taken from 
the same parent population. Categories which have P(null) < 0.1 for both tests are 
printed in bold. Values are based on all absorption features. 
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Table 8. Statistical comparisons of kinematic parameters: All galaxies 



Samples 


P(null, K-S) 


P(null, Kuiper) 


Av 


Seyfert Is vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


0.03 


0.09 


Seyfert Is vs Seyfert 2s 


<0.01 


<0.01 


AVm ax 


Seyfert Is vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


0.02 


0.02 


Seyfert Is vs Seyfert 2s 


<0.01 


<0.01 


Doppler Parameter 


Seyfert Is vs Starburst ULIRGs & LIRGs 


0.21 


0.03 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert Is vs Seyfert 2s 


0.02 


<0.01 



Note. — P(null) is the probability that the two listed distributions are taken 
from the same intrinsic distribution. Categories which have P(null) < 0.1 for 
both tests are printed in bold. Values are based on all absorption features, both 
inflowing and outflowing. Note that Seyfert 1 and Seyfert 2 galaxies include 
both IR-Faint Seyferts from the present study, as well as IR-Luminous Seyferts 
from RVS05c. Starburst data are taken from RVS05b. 
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Table 9. Statistical comparisons of kinematic parameters: Galaxies with outflows 



Samples 


P(null, K-S) 


P(null, Kuipcr) 


Av 


Seyfert Is vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


0.96 


0.84 


Seyfert Is vs Seyfert 2s 


<0.01 


<0.01 




Seyfert Is vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


0.26 


0.31 


Seyfert Is vs Seyfert 2s 


<0.01 


<0.01 


Doppler Parameter 


Seyfert Is vs Starburst ULIRGs & LIRGs 


<0.01 


<0.01 


Seyfert 2s vs Starburst ULIRGs & LIRGs 


0.24 


0.28 


Seyfert Is vs Seyfert 2s 


0.02 


0.01 



Note. — P(null) is the probability that the two listed distributions are taken 
from the same parent population. Categories which have P(null) < 0.1 for both 
tests are printed in bold. Values are based on outflowing absorption only. Note 
that Seyfert 1 and Seyfert 2 galaxies include both IR-Faint Seyferts from the 
present study, as well as IR-Luminous Seyferts from RVS05c. Starburst data 
are taken from RVS05b. 



Table 10. Statistical comparisons of kinematic parameters: Seyfert Is versus Seyfert 2s 

with inflows 



Samples 


P(null, K-S) 


P(null, Kuipcr) 


Av 


0.83 


0.87 




0.72 


0.89 


Doppler Parameter 


0.68 


0.99 



Note. — P(null) is the probability that the two listed 
distributions are taken from the same parent population. 
Values are based on objects with inflowing absorption 
only, and include data for IR-faint galaxies from this sur- 
vey and IR-bright galaxies from RVS05c. 
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Fig. 1. — Binned data showing the He I and Na I D complexes, plotted on an intensity versus 
wavelength scale, for four objects. In all cases, blue dashed lines correspond to the individual 
He I or Na I D line fits, with solid red lines showing the combined He I and Na I D fit for the 
object. Panels (a) and (b) show two Seyfert 1 galaxies, Mrk 506 and Mrk 6. Mrk 506 is fit 
with one He I and one Na I D velocity component. Mrk 6 is fit with two He I components 
and two Na I D velocity components. Panels (c) and (d) show two Seyfert 2 galaxies, Akn 
79 and NGC 3786. Akn 79 does not have a resolvable He I line; two Na I D components, one 
narrow and one broad, are used in the fit. NGC 3786 is fit with one He I component and two 
Na I D , one narrow and one broad. Note that the data for NGC 3786 are low signal-to-noise 
compared to most other objects, and thus the lines are not quite as well constrained. 
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Fig. 2. — Reduced x 2 values for free-floating and Ha-fixed He I parameter fits. See Section 
I4.ll for an explanation of these fits. The dashed black line corresponds to equal \ 2 values. 
The free fits generally give lower x 2 - The results from these fits were adopted in the present 
study. 



-38- 



> 
< 



400 



200 



-200 



-400 L 



-400 



-200 



1 1 1 1 

-♦Seyfert 1s 

- ■ Seyfert 2s 


/ 

/ 

/ 

✓ _ 


* 

— / 
/ 

/ 

, I I I I I I I 


I 1 1 1 I 1 1 1 



200 



Av free (km s ) 



> 

<3 



400 



000 



500 



o - 



-500 



000 



♦ Seyfert 1 s 
■ Seyfert 2s 



I 1 1 1 1 I 



500 t^i i i i I i i i i I i i i i I 
-1500-1000-500 



500 1000 



Av max: free (km s ) 



Fig. 3. — Left: Av for Ha- fixed He I fit versus Av for free-floating He I fit. Smaller points 
are objects with larger redshift uncertainty. Dashed black line through center corresponds to 
equal fixed and free values, and all points outside dashed black box have \Av\ > 50 km s _1 . 
Right: Same, but for Av max . The two types of fits give very similar results. 
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5840 5860 5880 5900 5920 

MA) 

Fig. 4. — Spectra of two Seyfert 1 galaxies, (a) Mrk 352 and (b) Mrk 6, in the Na I D region. 
These spectra were He I emission line subtracted, boxcar smoothed, and normalized to unity 
at 5910 A. The dotted line shows the 10/90 model (10%/90% stellar mass ratio for a young 
40 Myr instantaneous burst population compared to an older 10 Gyr population) and the 
dashed line shows the 1/99 model (l%/99% young to old stellar mass ratio), as outlined in 
Section 14.41 Strong stellar features are marked with dashed vertical lines. Mrk 352 shows 
absorption that is weaker than the stellar models, suggesting that the absorption is only 
stellar or atmospheric (the feature at ~ 5880 A was also determined to be atmospheric). 
In contrast, Mrk 6 shows very strong absorption, which is consistent with predominant 
interstellar origin. 




Av (Na I D P , km s" 1 ) 

Fig. 5. — Spectra of the He I emission + Na I D absorption region in infrared-faint Seyfert 2 
galaxies. These spectra are plotted on a velocity scale based on the systemic Na I D A5890 
velocity. The thin lines are smoothed original spectra, and thick lines are fits to the data. The 
vertical dotted lines show the locations of the Na I D AA5890, 5896 doublet and He I A5876 
emission line in the object's rest frame. 
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Fig. 5. — Continued 
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Fig. 6. — Plots showing the [N II] A6548 line for Av < and (properly scaled) [N II] A6583 line 
for Av > in 5 Seyfert 2 galaxies suspected of having blue emission-line asymmetry. The 
velocity scale is relative to systemic, and colors correspond to blue wing and red wing. The 
first four objects have noticeable blue-wing asymmetry, and the fifth (NGC 5728) shows a 
prominent blue bump in the emission line. The lower right panel (Mrk 622) is an example 
of an object which does not show an obvious asymmetry. 
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Fig. 7. — Spectra of the He I emission + Na I D absorption region in infrared-faint Seyfert 1 
galaxies. See Figure[5]for details. A region of horizontal lines is used to block out atmospheric 
absorption in the case of Mrk 352. 



45 




z = 0.0247 

I i i i i I i: i i i 



Akn_564 

i i i i i i i i i 



_' 1 1 1 1 1 1 1 1 




1 


1 1 1 1 1 1 1 1 1 1 '_ 


: z = 0.0293 

, 1 , , , , 1 ,: , , , 


i 


Mrk_10 : 

, , , i , , , , i , 




z = 0.0106 

i i i i i i 



1.5 



0.5 







I I 1 I I I 1 I , I I 


':' 1 1 1 1 1 1 1 1 1 '_ 


: z = 0.0335 

", , 1 , , , , 1 ,: , , , 


|Mrk_382 : 

, : , , , 1 , , , , 1 , ,~ 



i i | i i i i | i . i i i 


i . i i i | i i i i | i i 


z-0.0430 

", , 1 , , , , 1 ,: , , , 


iMrk_506 

, : , , , 1 , , , , 1 , 




z = 0.0169 

i .... i 



Mrk_359 

i i i i i i i 





1 1 1 1 1 1 1 1 1 1 1 '_ 


: z = 0.0374 

", , 1 , , , , 1 ,: , , , 


Mrk_477 : 

, , , , i , , , , i , ,~ 



i i | iii i | i.i 

a ■ ' 1 ^-i -• — ! 


1 1 1 1 | 1 1 1 1 | 1 1 


z-0.0274 

~, , 1 , , , , 1 ,: , , , 


Mrk_595 : 
, , , , i , , , , i , ," 



1.5 



■2000-1000 1000 2000 -2000-1000 1000 2000 



0.5 







Av (Na I D P , km s" 1 ) 



Fig. 7. — Continued 
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lo g( l fir / L © ) 

Fig. 8. — Histogram showing the fraction of Seyfert Is, Seyfert 2s, and starburst ULIRGs 
and LIRGs with outflows as a function of Lfir- The Seyfert categories include both IR- 
faint Seyferts from this study and IR-luminous Seyferts from RVS05c. Starburst ULIRGs 
& LIRGs are from RVS05b. The error bars are la, assuming a binomial distribution. Note 
the trend of increasing outflow detection rate with increasing Lfir for both starbursts and 
Seyferts. 
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Fig. 9. — Left: Histogram showing the percentage of velocity components as a function of 
Av for the IR-faint Seyferts from this study, the IR-luminous Seyferts from RVS05c, and 
the starburst ULIRGs & LIRGs from RVS05b. Objects at Av = —700 km s _1 actually have 
Av < —700 km s -1 ; the axis is truncated since most components have Av > —700 km s _1 . 
Solid black line at Av = highlights the division between inflows and outflows. Note the 
lack (excess) of outflowing (infalling) components in IR-faint Seyferts relative to the other 
objects. Right: Same as left, but the IR-faint Seyferts from this study and IR-luminous 
Seyferts from RVS05c have been combined together and divided into Seyfert types. Note 
the excess of infalling and high-velocity outflowing components among Seyfert 1 galaxies. 



< 

I 



CP- 

o 



1 1 1 1 1 1 1 1 1 

I ♦ Seyfert 1 s 


I I I I I I I I I 


- ■ Seyfert 2s 




: •Starburst ULIRGs 


& LIRGs I 


7X Starburst Dwarfs 






♦ - 

X - 


XX x 









5 



1 .0 



S 3 

E 



< 

I 



CP 

o 



2 - 



T 



♦ Seyfert 1 s 
■ Seyfert 2s 

• Starburst ULIRGs &c LIRGs 
- Starburst Dwarfs 



1 .5 

log ( v c 



2.0 



2.5 
1 



3.0 



/ km s ) 



XX X 



9 10 11 12 

log ( L FIR / L ) 



1 3 



Fig. 10. — Left: Maximum outflow velocity versus circular velocity. The Seyfert 1 and Seyfert 
2 categories include data from this study and from RVS05c (both IR-faint and IR-luminous 
Seyferts). The starburst ULIRG & LIRG category is composed of data from RVS05b. The 
data on the starburst dwarfs are from Schwartz & Martin (2004). Right: Same as left, but 
for maximum outflow velocity versus Lfir- The outflow kinematics in Seyfert 2s (Is) are 
similar to (different from) those of the starburst ULIRGs and LIRGs. 
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Fig. 11. — Left: Mass outflow rate versus circular velocity. See Figure fTUl for descriptions of 
the various categories. Right: Same as left, but for mass outflow rate versus L FIR . 
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Fig. 12. — Left: Av max versus circular velocity for inflows detected in IR-faint Seyfert 
galaxies. Smaller points are those with uncertain redshifts. Right: Same as left, but for 
Av max vs. Lfir- No obvious trend is seen. 



